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I. INTRODUCTION 



The expansion of matrix- valued field theories is probably the most important 
non-perturbative and non-numerical theoretical tool presently available in the study of such 
models as non-Abelian gauge theories and two-dimensional quantum gravity. A resolution of 
the above-mentioned models in the large- limit would be the starting point for many ana- 
lytical developments. In particular when a lattice formulation is involved one must consider 
different possibilities in the search for the continuum limit; for the case of asymptotically 
free theories one must explore the limit of vanishing coupling g = (trivial fixed point) while 
keeping a physical mass scale fixed, while in the case of quantum gravity one must search 
for a nontrivial fixed point Qc and reach the limit with a specific power-law dependence on 
A^ of g — gc, which is known as "double scaling limit" [0,0,0. Therefore it is useful to achieve 
a full knowledge of the coupling dependence of such models, from extreme weak coupling to 
strong coupling, in order to explore those regions that may turn out to be physically most 
interesting. 

As a matter of fact, notwithstanding many recent efforts toward an understanding 
of the possible properties of large- A^ solutions to nontrivial quantum field theories, our 
present analytical knowledge is limited to a small number of few-matrix systems. This 
number is even smaller if we restrict our attention to the case of unitary matrix fields, which 
is especially relevant to the problem of lattice QCD. To the best of our knowledge, the 
only solved examples are Gross- Witten's single-link problem and its generalizations, the 
external field problem and L = 3,4 chiral chains [^|[. 

We stress that extending the number of solved few-matrix systems is not at all a 
pointless exercise. Indeed apart from purely theoretical informations that might be achieved, 
not only does every few-matrix system have a reinterpretation, via the double scaling limit, as 
some different kind of matter coupled to 2-dimensional quantum gravity, but also every few- 
matrix system involving unitary matrices can be reinterpreted as the generating functional 
for a class of integrals over unitary groups, and these integrals in turn are the essential 
missing ingredient in the context of a complete algorithmization of the strong coupling 
expansion of many interesting models 0. 

Following Ref. we may introduce the notion of a "superskeleton" , that is a graph 
whose vertices are joined by at most one link (simple graph). As has been shown, knowledge 
of all the group integrals involved in the strong coupling expansion of a lattice model with 
nearest-neighbor interactions defined on such a graph provides sufficient information for the 
algorithmic reconstruction of the strong coupling series for a model enjoying the same global 
symmetry and defined on an arbitrary lattice. 

These were basic motivations for us to begin the study of the class of lattice chiral 



models which we termed "simplicial chiral models" |10]. In particular we focused on principal 
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chiral models, with a global U{N) x U{N) symmetry, defined on a ci — 1 dimensional simplex 
formed by connecting d vertices by {d — l){d — 2)/2 links, and explored specifically the 
large- limit of such models, whose relevance we have just been discussing. 



Our fundamental result is the reduction of the above problem to that of solving a 
single inhomogeneous integral equation for the eigenvalue distribution of a single Hermitian 
semi-positive definite matrix. Although we could not find a closed form solution to this 
equation for arbitrary d, we are able to solve it in several interesting special cases and we set 
up a systematic numerical approach to the solutions which led us to a conjecture about the 
location of the critical surface as a function of d. We have also studied in detail the related 
topic of chiral chains, their strong coupling expansion and critical behavior. As a result of 
these analyses, we are confident that the critical surface is defined by /3c = ^/d for all d. 



Moreover, by treating d as a continuous parameter, there are two distinct regions. For 
small d, < d < 4, the models exhibit the third order Gross- Witten transition. Indeed 
for d = 1, 2, 3 they coincide exactly with the chiral chains studied earlier by Brower, Rossi 
and Tan |^. In this region, the criticality is related to that of 0{n) spin models on random 
surfaces, as discussed by Gaudin-Kostov [|T3|. For d > 4, however, there is a first order 
transition ending at the "upper critical" dimensions d = 4, which we scrutinize in some 
detail. 



This paper is organized as follows: 



In Section || we set up our formalism for simplicial chiral model and derive the large- A^ 
effective action and a representation for the internal energy. We begin in Section WTB 



by illustrating the formalism for the simpler case of vector 0{N) spin models on a 
simplicial lattice, deriving the closed form large-A^ solution for arbitrary values of d 
and studying its properties. Sec |11 C| gives the large- A^ effective action for the simplicial 



chiral model and Section |II D| , the saddle-point equation (large- A^ Schwinger-Dyson 
equation) for the eigenvalue distribution, discussing its features and converting it into 
a standard inhomogeneous Fredholm equation of the second kind. 

In Section |T| we analyze the solvable examples. We begin with integer values of < 4 
(which correspond to L < 4 chiral chains), followed in Section |III B| by a detailed 
discussion of the d = 4 case, including features of the weak and strong coupling expan- 
sions and the asymptotic expansion around the critical point jSc = 1/4. We present 
in Section [III C| the exact result for the limit d ^ oo and develop in Section [III D| a 



treatment based on the 1/d expansion. This large-d analysis, which works best outside 
the critical region: j3 ^ jSc, has provided us with the first numerical confirmation for 
the conjecture that jSc = 1/d, and is discussed in greater details in Appendix D. 

In Section |^ we solve the models at criticality for arbitrary < < 4 and sketch 
the peculiar features of the critical behavior when d > 4. Also in Section [IV B[ we 
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present numerical methods based on the Gaussian integration techniques for d > 4 on 
the critical surface (3c — 1/d. 

• Some technical extensions are included in other appendices. Appendix A is devoted to 
a discussion of the double- scaling limit of critical chiral chains L < 4 and Appendix B 
extends the discussion of chiral chains to L > 4 by an analysis of strong coupling 
expansion. Whereas Appendix C is devoted to details of the weak and strong coupling 
expansions and series analysis for the d — A simplicial chiral model. 



II. SIMPLICIAL CHIRAL MODELS 



A d — 1 dimensional simplex is formed by connecting in a fully symmetric way d 
vertices by {d — l){d — 2)/2 links. Let us assign a U(N) matrix to each vertex. The partition 
function for principal chiral models on a simplicial lattice is obtained by integrating over 
unitary N x N matrices with a normalized invariant Haar measure: 



= [lldUi exp In f3 ^ Tr [UiU] + UjU[ 

^ 1=1 I i>j=l 



(1) 



Thus the d-matrix simplicial model has an underlying permutation symmetry instead of the 
cyclic symmetry of the d-matrix chiral chains. For d— 1,2 and 3 these two symmetries and 
the associated models are equivalent. We shall explore the U{N) x U{N) symmetry of the 
system and, in particular, study its critical behavior in the large- A?" limit. 



For this purpose, it is sufficient to study the bulk "thermodynamic" properties, e.g., 
the free energy density, internal energy, and specific heat, which are respectively given by 

FNiP^d) = ^In^d, 

C^(f),d) = fi^^S. , (2) 



We focus in this paper on computing the free energy and determine the critical point (3^ for 
all values of the parameter d in the large A^-limit. We find that there is a sequence of critical 

theories with I3c = 1/d, which exhibit a third order Gross- Witten singularities for < 4 and 
a first order transition for d > 4. Special attention will be give to the marginal dimension at 
d = 4. Scaling exponents and finite size (or double scaling) properties will be presented in 
some special cases, but a thorough investigation for all d is beyond the scope of this paper. 
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A. External Field Method 



Because of the the permutation symmetry of the vertices, the simphcial chiral models 
can be reformulated in terms of a Lagrange multiplier field which decouples the original 
degrees of freedom. The resulting effective theory is very reminiscent of the mean field 
approximation to standard lattice models, but in contrast with mean field this reformulation 
is exact. 



We therefore replace the direct interaction between unitary matrices with the coupling 
to an auxiliary field, which in this case is a complex N x N matrix. A, introduced in the 
following representation of the identity. 



/ dA exp {-Nf3 Tr [{A - Eti Ui) (^t - Ef=i U^)] } 



JdAexp{-N(3TT [AA^} 



(3) 



Again by exchanging the order of integrations and representing the partition function in the 
form Zd = Zd/Zo, we may obtain 



Zd = dA II dUi exp \ N/3Tt 



1=1 



AA^' 



+ 



+ 



d 



(4) 



By performing the single-link external field integral, we may introduce the auxiliary function, 
F^(SSt) = ^In J rff/expj^Tr [S[/t + [/St]| , (5) 



and re-express Zd, up to an irrelevant multiplicative factor, in the form 



dB exp 



-^TiBB^ + N^dF(BB^) - N^(3d 
A/3 



(6) 



where B replaces 2(3A. 

A first crucial point in our analysis is the observation that the integrand in Eq. (^) is a 
function of the eigenvalues Xi of the Hermitian semi-positive definite matrix BB'^ . Moreover 
Morris |T2| has shown that, when integrating over complex matrices, a proper parameter- 



ization may offer in specific cases, like ours, the possibility of performing the "angular" 
integrations exactly and reducing the problem to that of integrating over the N variables Xj. 
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Referring to Morris' paper for a proof of the angular integration, we apply it to our Eq. 
thus obtaining, again up to irrelevant numerical factors, 

_ poo d. r 

Zd = l[dxiY[{xi-Xjfexp-—J2xi + N^dF{xi)-N^pd 



(7) 



The second crucial observation concerns the function F{xi). It is known exactly for all 
U{N) groups [0], while integral representations exist for SU{N) groups [0, and it takes 
on a relatively simple form in the large- limit. Before proceeding further, we shall first 
provide an even simpler illustrative example whose large- solution can be obtained fairly 
straightforwardly. 



B. Pedagogical example: simplicial spin models. 



Consider instead of our simplicial chiral models, an example of an 0{N) symmetric 
nonlinear model defined on a simplex. The same basic methods used for the chiral models 
are easily illustrated in this much simpler context. 

The partition function is obtained by integrating over the — 1 independent compo- 
nents of d vectors, 



Zd = Jll dsi 5(s2 - 1) exp 

i=l 



i>j=i 



(8) 



The effective field is a single unconstrained A^-component vector, which again can be 
introduced as Lagrange multipher field via an identity. 
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/rfzexp-^(z-EtiS. 
/c/zexp-^z2 



(9) 



Upon substituting this identity into Eq. (H) and inverting the order of the integrations, we 
may then represent the partition function by = Z^/Zq, where 



i=l 



Zd = dvY[ dSi 6{s^ — 1) exp 



.fv^...v.(|:.)-f. 



(10) 



It is now possible to perform the decoupled constrained integrations. To this end we may 
define the auxiliary function 
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In y (is S{s^ — 1) exp A^s ■ z 



(11) 



where = z ■ z. This function is known exphcitly for all values of N and admits a large N 



limit 1 1 



liur(f)/„_.(iv. 



1 



- 1 - In ( -Vl + 4^2 + - 



(12) 



As a consequence for large N we have the following representation of Zd, 



dzexp < — - 



- - (^Vl + 4^2 - 1 - In (^-71^412 + -)) + 



(13) 



The large- A^ value of the integral in Eq. (|13]) may be obtained by a saddle-point estimate. 
After some simple manipulations, the saddle-point equation may be reduced to 



i-5^|(i-v/rT4i^). 



(14) 



The solution of this equation is 



l-^ + 4 + ^V(l-/5^)' + 4/3 



(15) 



By taking the logarithmic derivative of the partition function with respect to (3, we obtain an 
expression for the internal energy (per unit link) Ud of simplicial spin models in the large- A^ 
limit 



d{d-l) _ I fz^ , 1^ _ 



4/3-1 



(16) 



We may check many special cases of this result, and in particular we may notice that the 
r.h.s. of Eq. (|TB]) is zero when d = 1, while when d = 2 



U2 = — 

' 2(3 



1 + 4/32-1 



(17) 
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consistent with the single-hnk model result. 



Finally let us notice that in the large-d limit, as a trivial consequence of the structure 
of the model, the solution we found coincides with the mean field solution, which is exact 
in this limit. It is worth observing that, while Eq. ([T6|) is formally correct for all values of 
/?, in order to recover the standard strong and weak coupling expansions of the solution we 
must separately consider the two different regimes [3 < [3c and (3 > Pc, where for all d we 
obtain (3c = 1/(1. 



C. Large- limit 



Returning to the simplicial chiral models, we are again interested in the large- limit. 
For the free energy function F{xi) resulting from a one-link integral over a U{N) matrix, 
the limiting form can be extracted by solving the Schwinger-Dyson equations and written 
in a simple closed form f^J^, 

s 1 r-, 1 f Vr + ^i + Vr + Xj \ r 3 



We must distinguish two different phases, a weak coupling regime where r = and 

—y^ < 1, (19) 



and a strong coupling regime where r is dynamically determined by the condition, 

- y , ^ = 1 • (20) 



It is important for future developments to observe that Eq. (^) also leads to the condition 



dr 



. (21) 



It is completely legitimate to apply the above results to a saddle-point evaluation of the 
large- limit of the integral appearing in Eq. (0). To this end we may define an effective 
action 
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Sd = 775 E ^« - N'dF{x,) - ^ ln(x, -x,) + N^(3d , (22) 



and derive a saddle-point equation 



Very simple manipulations, including the use of Eq. (pT]), lead to a reformulation of Eq. (p3D, 
which can be turned into the relationship 

Vr + Xj _ . _ J_ (4 - d)^/r + Xi + dy/r + xj , , 

2/3 " iV^ x,-x, ■ ^ > 



This equation, supplemented with the condition > and with the constraint r = (weak 
coupling) or Eq. (|20|) (strong coupling) is the fundamental saddle-point equation of principal 
chiral models on a simplicial lattice. It is the starting point of most of the developments 
presented in the following sections. 

We recall that, once Eq. (^) is solved, knowledge of the saddle-point value of the 
eigenvalues Xi allows the large- evaluation of Zd via the relationship, 

Zd — > exp[-S'd(xi)] , (25) 



and we can also extract the internal energy per unit link by taking a logarithmic derivative 
of Zfi with respect to (3 which leaves us with the relationship, 

d{d-l)U = ±-^J2^^-d-^■ (26) 



D. Saddle-point equation. 



In order to study Eq. (^J) we shall start by applying well-established techniques, and 
in particular by introducing an eigenvalue density function. It is however convenient first to 
introduce a new variable Zi whose formal definition is 

Zi = Vr + Xi. (27) 
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subject to the condition < -Jr < Zi. We may assume that the eigenvalue variable Xi lies 
in a single interval [xa,^;,], < < Xf,. In terms of the new variable one has 2je[a,6] 
where a = ^Jr + x^, h = + x^ and 

< < a < 6. (28) 



For weak coupling, r = 0, and we expect in general a = ^Jx^ > 0. For strong coupling, one 
expects Xa = so that a = -y/r 7^ 0. 

We shall be interested in the weak-strong transition as one varies d and f3. In a 
third-order transition, typical of large-iV transition previously studied, Oc = v^c ~ 0- -'■^ ^ 
first-order transition, which we will encounter for d > 4, = \fr^ 7^ when approached 
from the strong coupling regime. 

Denoting the large- eigenvalue density by p{z)\ it vanishes outside the interval [a, 6]. 
We may now turn Eq. ( P^ into the following integral equation, 

— -d= J^'dz'piz') ( ^— - . (29) 

2(3 Ta ' \z-z' z + z'J ^ ^ 



The function p{z), and therefore also the extremes a and b of the integration region, are 
thus determined dynamically. In particular the normalization condition, 

rb 

dz'p{z') = 1 , (30) 



must be satisfied. 



In addition to the positivity requirement, p{z) > over the interval [a,b], the desired 
solution to Eq. (|29| ) must also satisfy either the weak coupling inequahty, Eq. ([I9|) or the 
strong coupling constraint, Eq. (pOl). In the large- limit, Eq. (^) becomes 

'dz'^ < 1, (31) 



whereas Eq. (pO]) becomes 



'd.'<!!^ = l. (32) 



The determination of the transition point, /?c; and of the critical behavior around this value 
is one of the interesting physical problems concerning this model. 
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Eq. (^) has a somewhat unconventional form when compared to other integral equa- 
tions, because of the special structure of its kernel. We may however perform a few manipu- 
lations in order to obtain a more familiar relationship. Our starting point is the introduction 
of an analytic function oi z, by the definition 

fiz) = f^dz'. (33) 



By construction, the analyticity domain of f{z) is the complex z plane with the exception 
of a cut on the positive real axis in the interval [a, h]. The discontinuity on the cut may be 
parameterized by writing 

f{z ± ie) = R{z) T iTrp{z) , (34) 



when z G [a,b] and it is easy to recognize that 

^ , , z d d — 2f^. p(z') z d d — 2^, . , ^. 



It follows that R{z) is itself an analytic function of z, with a cut on the negative real axis 
in the interval [—6, —a]. 

Let us now notice that the normalization condition implies 

IzHoo Z 



As a consequence in the same limit we obtain 

R{z) ^ A _ ^ + IZI^ (37) 



and in turn 



1 / z d d — ' 
V 4/3 ~ 2^ ^7 



O 



(38) 



This equation can in principle be used in order to determine relationships between the 
constants a and b in place of the normalization condition. 
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We must now distinguish between weak and strong coupling regimes. In both cases, 
by exploiting analyticity properties of the function f{z) and defining appropriate auxiliary 
functions, it is relatively easy to reduce Eq. (p9D to the following forms 



\jip-z){z-a) 



TT 



1 

4^ 



d-2 



dy 



piy) 



- y + z^{h + y){y + a) 



for /3>/3e, (39) 



Tt \ Z 



1 

4^ 



d-2 



dy y + ap{y) 



y+z^y+b y 



for P<Pc 



(40) 



The values of a and b as functions of /3 are determined by enforcing the asymptotic condition 
(^). Eqs. (|39| - ^0|) are inhomogeneous Fredholm equations of the second kind. It is therefore 
in principle possible to apply standard methods of (approximate) resolution by expressing 
the kernels in terms of appropriate orthonormal sets of eigenfunctions. 



The weak and strong coupling constraints, Eq. (|I^ and Eq. (^Oj), can be expressed in 
terms of the analytic function f{z) as /(O) > —1 and /(O) = —1 respectively. Alternatively, 
writing f[z) = R{z)~mp{z) and analytically continue this expression outside of the interval 
[a, 6], Eq. ( pj]) and Eq. (pO|), can also be expressed as 



-ip(0) > 



(41) 



and 



-ip(0) = 



(42) 



respectively. Note that Eq. (|40|) is parameterized so that the strong coupling constraint 
(^) is automatically satisfied. The transition point (3c can be determined approaching from 
the weak coupling regime by enforcing the equality p(0) = 0. We shall return to a general 
discussion of this criticality in Section IV. 



A final comment concerns the explicit evaluation of Zd- Instead of directly substituting 
p{z) in the expression of the partition function, it is convenient to apply Eq. (|26|) in the 
form 



d{d-l)U 



1 



4/32 



b I 

dz' p{z){z''^ — r) — d — — 

P 



(43) 



and perform an integration with respect to (3 to recover the free energy. 
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III. EXACT LARGE N SOLUTIONS 



Here we present the solutions at = cxd as a function of d. They can be broken into 
three classes for d < 4, d = 4 and d > 4 respectively. For d < 4, they are equivalent to chiral 
chain models with L < 4 studied earlier all of which exhibit the third order Gross- Witten 
transition at (3c- For d > 4 there is a first order transition, which ends exactly at d = 4, 
consequently the end point at = 4 is of special interest. 



A. Solutions for d < 4 



As we mentioned in the introduction, when d is integer and less than 4 simplicial 
chiral models are only reformulation of trivial or already solved models. It is however quite 
instructive to consider even these examples in our new language. Let us begin with the only 
apparently trivial case d = 0. Obviously Zq = 1, however Zq is nontrivial and we need to 
know its value in order to compute Zd- As a matter of fact Eq. already implies that, up 
to a constant 

Zo oc expNHn(3 . (44) 



We would like, as a consistency check, to derive this result from the saddle point equation. 
A straightforward manipulation of Eq. (^) leads to 



Ja 



— z' 



This equation is solved by 



with the only constraint IP' — o? = 16/3. However by keeping in mind that only the combi- 
nation X = z^ — a? is physically meaningful because of Eq. (^), we recognize that 



piz)dz = — , = — -^Jl6p-x , (47 

'^^ ^ 877/3 _ ^2 87r/3v^V ^ ' V ^ 



and the physical solution is unique and leads by a trivial integration, to Eq. (pj). 
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Because of our definitions Z\ = 1, Zi = Zq. The eigenvalue distribution p{z) however 
is the generating function for the moments of the linear combination of a complex and 
a unitary matrix, and these moments can be highly nontrivial, even if the complex matrix 
itself has a Gaussian probability distribution, as a consequence of the averaging over unitary 
matrices. As a matter of fact we were not able to solve explicitly the saddle-point equation 
associated to the d = 1 models, even though the solution probably has reasonably simple 
mathematical properties. 

Let us now turn to the d = 2 case. As a straightforward application of Eqs. (39|-40), 
we immediately find both the weak and strong coupling solutions, 

p^{z) = J^^8f3-{z-A(3y for /5 > ^ , (48) 



Psiz) = + for /3 < 1 . (49) 

^ ^ ^ Att(3 \ z -1 + 2(3 ^ - 2 ^ ' 



For the strong coupling region, r = (1 — 2/3)^. 

It is easy to recognize that the d = 2 model corresponds to the Gross- Witten single-link 
problem, which in turn is equivalent to large- QCD2 with Wilson action on the lattice . 
The properties of this model are well known, and in particular it is known that (3c = 1/2, 
consistent with Eqs. (|48| - ^9D . Another consistency check is easily made by applying Eq. (^31) 



and verifying that the known expressions for U2 are reproduced. 

Finally let us comment about the d = 3 case. This model in its original formulation 
is completely equivalent to the three-link chiral chain studied in Refs. We therefore 

know that it must possess a third-order phase transition at the critical value (3c = 1/3. 
However, as we already observed, our reformulation leads to exploring quite different classes 
of correlation functions and there is no obvious relationship between old and new results 
apart from bulk thermodynamical properties. Again we have no analytical solution for the 
d = 3 model equation, whose known properties stand as a benchmark for future attempts. 



B. Solution at d 



Turning to c? = 4 leads us to a new situation, where we are no longer guided by known 
results, since the 3-dimensional simplex (tetrahedron) is distinct from the solved four-link 
chain. Actually it would be instructive and convenient to embed both models in a more 
general case interpolating between them and including many more interesting situations. 
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We are studying the most general four-site system with bihnear interactions of four unitary 
matrices, which turns out to be reducible to an interacting two-complex matrix system. A 
separate paper will be devoted to a discussion of this system. Here we only discuss the 



solutions of the saddle point equation obtained from Eq. (pQ]) in the d = 4 case, 



8/3 



In order to solve this equation, let us separately consider the weak and the strong 
coupling regimes, while changing variables for convenience to x = z"^ and defining the dis- 
tribution p{x) by p{z)dz = p{x)dx. The special structure of Eq. (|50|) makes it convenient to 
follow a special procedure not directly related to Eqs. (|39|-^0D derived for the general case. 



In the weak coupling phase, we define the functions 



U{x) = I ^ dx , (51) 



and 



g^x) = , ^-^""^ , (52) 

{x — lP'){x — a?) 



subject to the normalization constraint 



dx'p^{x') = 1. (53) 



The functions f{x) and g{x) are real analytic, with a cut along the interval on the 

real axis. On this interval the relationship, 

lmg^{x±ze) = T , ^"^"^""^ (54) 

(x — a'^){lP' — x) 



holds, while analyticity and Eq. (pOj) imply 



8/3 1 ^ r , Vy 



9wix) = I - 7-— 7 / dy — , . (55) 

^{x-h^){x-a?) 8npJo ^ + 2/ ^^(2/ + 62)(^ + a^) 
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However Eqs. (pl|-|5^) imply that 



lmf„{x + ie) 



dy [y 



tta/x 8n^j3 Jo x + y\ {h^ + y){y + a 



(6^ — x){x — a^) 



In order to determine o? and 6^ we may use Eq. (|53|) and the observation that in 
X plane when |x| ^ oo 



x" 



As a consequence one obtains that 

j) dx g^{x) = 2j ^ dx'lmg^{x') = 



around [a , 6 ], that is 



fe2 



dx' 



{b^ - x'){x' - a^) 



^71(3 



In strong coupling we adopt a similar strategy by defining 

r^' dx' ps{x') 



fs{x) 



a2 A/a/ X — x' 



and 



9s{x) = /s(a;)i 



/X — 
X — h"^ 



with the constraint 



I? 



Ps{x')dx' = 1 



and the boundary condition 



9s[x) 



1 

\x\^oo X 
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We then find 



1 X — a? 



9s[x) 



dy 1 jy + 



(64) 



and 



62 



dy x{b'^ — x){y + a^) 
Sn'^lS X + y\ yilP' + y){x — a^) 



(65) 



The boundary condition leads to the relationship 



— J dx' Im gs{x' + ie) = vr 



(66) 



where 



lmgs{x + ie) 



1 X — a'^ 



8/3 a/x ^ IP' — X 



(67) 



All the integrals appearing in Eqs. (|56|), (|59D, (|65|) and (|66|) are elliptic integrals. It 
is therefore possible to re-express both the weak and the strong coupling results in terms of 
known functions. In particular it is convenient to re-express everything in terms of "natural" 
rescaled variables, by defining 

(68) 



(69) 





and setting p{C)d( = p{z)dz. It is not too difficult to eliminate completely the parameters 
a and b in favor of k by making use of Eqs. (^) and (|6^) respectively. As a consequence we 
obtain the weak coupling expression 



Pw(C) 



E{ky 



(70) 
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and the strong coupling counterpart 

Ps(C) 



[E{k) - (1 - k'^)K{k)Y 



(71) 



where K,E,Il are the elhptic integrals of the first, second and third kind respectively, and 
the domain of ( is the interval [0, k], < k < 1. 

Obviously, in order for the problem to be completely solved, one must try expressing 
A; as a function of f3. This is achieved in principle by enforcing the normalization condition, 
which takes the form 

f-p{OdC = l. (72) 
Jo 



By symbolically writing 



p(C) = mC,k), (73) 



in agreement with Eqs. (^) and (0), it is actually possible to express all results as functions 
of k by the relationship 

^ ^ J^kWJ) ■ ^^^^ 



In practice this form of our results is sufficient for both numerical evaluation and 
asymptotic expansions, not to mention the possibility of exploring the region around the 
criticality. Criticality is characterized by the limit A; — 1, where simple mathematical 
properties of elliptic integrals allow us to show that both weak and strong coupling results 
lead to /?c = 1/4 and 

Pc(C)=Cln[^. (75) 



In order to obtain the usual weak and strong coupling expansion of physical quantities, 
like the internal energy, as power series ml/ (3 and (3 respectively, one must consider in turn 
the A; — >■ limit and the expansion in powers of k. Obviously the different structure of p{Q 
in the two phases will lead to different expressions. In particular we have the asymptotic 
behaviors 



18 



k^o 16 32 



(76) 
(77) 



and it is conceptually straightforward to obtain power series expansions in the powers of 
k for such quantities as the internal energy and to convert them into standard weak and 
strong coupling series. A few details will be discussed in Appendix C. 

The expansion around the critical point /^c = 1/4, fc^ = 1, is slightly subtler because 
the expansion of elliptic integrals around A; = 1 is asymptotic. However by exploiting a few 
known or previously derived results, we have managed to obtain the following relationships, 
holding in weak coupling near the criticality: 



Pw(C) ^ 4/3 



l-k'^{ In 



Oik 



I2\ 



Cln 



1 + C 
1-C 



k'^ ( In 



1\ 



k' 27 1-C 



+ 0{k 



I2\ 



(78) 



where k''^ = 1 — A;^ ^ 0. From Eq. (|7^ we then obtain 

[42121 

4/3 ^ l + k'"^ In — In h - In h - 

^ L A;' A;' 2 A;' 2 



(79) 



and as a consequence 



k'^ [ 4 2 1 2 1 



(80) 



therefore k''^ ^ (3 — /3c apart from logarithms. By properly applying Eq. (^) we may also 
extract the result 



7r^-6 

k'^O 9 



+ Oik'Hn^ k'^ 



(81) 



and by simple manipulations, from the specific heat relationship 



(82) 



we may obtain near the criticality 
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7r2 + 3 



vr 



O 



k'^o 36 121n(4/A;') Vln^A;' 



(83) 



A similar analysis can be performed in strong coupling near criticality, 



Ps(C) ^ 4/? 



4 1 



Cln 



1 + C 
1-C 



A;'2 ( In _1 + 1 



4 , 1\ C 



k' 2/ 1 - C 



0(A; 



/2\ 



(84) 



where again k 
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1 — A;^ ^ 0, and 



J2 



4 2 
In — In — 



2 A;' 2 



(85) 



We then find 



TT 



TT 



fc'-o 36 121n(4/A;') Vln^fc' 



(86) 



The strong and weak coupling expressions of C near the criticality show that the critical 
behavior around [3c = \ corresponds to a limiting case of a third order phase transition with 
critical exponent of the specific heat 



a = 0" 



(87) 



near the boundary with weak second order critical behavior. Notice that in terms of double 
scaling limit a = 0^ would correspond to a central charge c = 1. 

For the interested readers we mention that in the derivation of Eqs. ([78| ) and we 
made use of the following formula (which appeared with some misprints in Ref. |ir 



^ 4 C, l-C 
In — + - In ^ I + 



1 - CM k' 2 1 + C/ 4(1 - C 



2^2 



(88) 



for the asymptotic expansion of the elliptic integral of the third kind 11 (C^, k) in the region 
k' = v^r^F ^ 0. 
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C. The d = oo solution 



While at present we are not aware of any general method to get an analytic solution of 
the saddle-point equation (^) for arbitrary d, the d —>■ oo limit provides another interesting 
instance in which the equation is solvable. 

It is easy to show that for larger and larger values of d the distribution p{z) becomes 
narrower and narrower, with a width decreasing like d"^/^ and a peak value z which can 
easily be determined by replacing in Eq. (^Of ) 

p{z) S{z - z) , (89) 



and obtaining the large- A^, large-rf equation 



z , d — 2 , . 

— -d= . 90) 

2P 2z ^ ^ 



A consistent solution is obtained by assuming the limit to be taken at a fixed value of f3d, 
in which case 



with the obvious restriction jSd > 1 (weak coupling phase). When l3d < 1 one must recognize 
that the saddle-point condition, when correctly applied to the original expression for the 
effective action Eq. (^) in the large-ci limit, unambiguously leads to the prediction z ^ 1, 
r — > 1 (strong coupling phase). 



The most interesting features of this result are: 



• The large-d prediction for the location of the critical point, (3c — ^/d, amazingly 
enough, seems to be satisfied for all values of d. 

• The complete equivalence with the mean field solution of infinite volume principal 
chiral models on a D-dimensional hypercubic lattice such that D = d/2 |jT^, where 
we may observe that this last relationship enforces the constraint that corresponding 
models have the same coordination number. 



It is easy to compute the large-c? expression for the internal energy in the weak coupling, 
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At the criticality. the weak couphng value of [/ is ^, while the strong coupling value is 
U = 0. Therefore the large- A^, large-d prediction for the nature of the criticality is that of 
a first-order phase transition. It is however important to notice that the large-ci prediction 
for the specific heat in the weak-coupling phase, 



dC 



1 



^ I3d 



+ 1 



(93) 



shows a divergence at the phase transition, with no indication for the existence of a 
met ast able phase. 

It is interesting to compare the specific heat behavior for d = 2,3, 4, oo. In Fig. |I| we 
plotted dC versus (3 = dp. 



D. The 1/d expansion 



The large-c? result may also be the starting point for a systematic 1/d expansion of 
Eq. (^), and for a numerical approximation scheme which turns out to be quite efficient at 
least in the weak coupling domain away from criticality. The essential ingredient for both 
these developments is the observation that, substituting the definition of f{z), Eq. (|3^), into 
Eq. (p9|) we obtain the functional equation. 



2/3 



-d 



2Re f{z) + {d- 2) f{-z) 



(94) 



subject to the following constraints: (a) Eq. ( pi] ) is satisfied in the interval [a,b], with a 
and b dynamically determined; (b) f{z) is real analytic outside the interval [a,b]; (c) the 
asymptotic behavior of f{z) when \z\ — > cxd is f{z) — > 1/z. 



Let us now define 



C^z-~z, 0(0^/(5 + 



(95) 



where z is a constant whose value lies in the interval [a, b] and will be dynamically generated. 
Substituting the definition (031) into Eq. (Pl), we obtain 
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^-d + ^ = 2Re 0(C) + {d- 2)<P{-2~z - Q . (96) 



Postponing the discussion of the numerical approximation scheme, let us illustrate here the 
procedure for a systematic 1/d expansion of Eq. (pG]). We introduce the following Ansatz 
for the function 0(C), 



m = dQiC) [I - ^1 - ^ - :^ ) + ^(C) , (97) 



dC dc 



and assume the functions Q{C), -R(C) to be real analytic in the interval between the roots of 
the polynomial 

dC' - ciC - C2 , (98) 



including the point C = 0. Moreover we require Q{0) = R{0) = and the boundary condition 

m , ^ . (99) 



These requirements fix the constants z, Ci and C2 dynamically, which in turn determine a 
and b. Finally we assume all functions and constants to be expandable in 1/d, with non 
vanishing leading order. 

As an illustration let us find the solution to first nontrivial order. All leading order 
quantities will be independent of d and labeled by a subscript 0. After an expansion of 
Eq. ( pTl) in powers of l/d we obtain 

0(0 ^ ^Qo(C) (f + ^) + ^o(C) + oQ . (100) 



By imposing the asymptotic boundary conditions we have the stricter condition (forced by 
analyticity of Qo, Rq), 



iQoiC) [j + + RoiO = ^. (101) 



As a consequence we may also predict the O (^2) asymptotic behavior. 
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(102) 



Now by substituting the above results into Eq. (^) we obtain, after expansion in powers of 
1/rf, 



(103) 



where we always assume the large-d limit to be taken while keeping fid finite. Substituting 
the condition Qo(0) = into Eq. ( |103| ) we may solve it in the form, 

1 11 

(104) 



2f3d zo 2~zl ' 



Qo(C) 



25o 



+ 



c 



45o/ 42o25o + C 



(105) 



while the implementation of Eq. ( |101|) fixes 

1 1 

ci -- 



C2 



2^0 f 1 _ ^ 

V 420 

4zo 



(1 



3 

420 



2 ' 



(106) 



and in conclusion we may also write 

^o(C) = 
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850^1 -it)' 2'=° + « 



(107) 



Let us notice that the leading order is completely determined in terms of the parameter So, 
which in turn is fixed through Eq. (|104|) to take the value 




5o = z = /5ci 1 + Jl-— . (108) 



Hence we recognize that z is nothing but a generalization of the mean field parameter which 
(roughly speaking) describes the center of the eigenvalue distribution, while the width of the 
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distribution itself is O (l/Vd) as one may easily see by studying the roots of the polynomial 
under the square root sign. 

The eigenvalue distribution itself may be recovered (order by order in 1/d) by the 
relationship 



dQ{z - z) jC2 ci 



p{z) = ^ + ^{z -~z)-iz- ~zY , (109) 

vr U — z) \ a a 



and it is not too difficult to check that the large-d limit of Eq. ( |109|) may be taken and the 
result is 

d 2 



lim 



C2 
d 



- u 



5{z- zq) 



(110) 



as expected. We have worked out higher orders of the 1/d expansion. 



The Ansatz (pTl) can also be used as a starting point for numerical approximations, 
and this will be described in Appendix D. For instance, with the precision of about 1%, one 
can quickly verify numerically the conjecture that j3c = Xjd. 



IV. SIMPLICIAL MODELS AT CRITICALITY 



We would like to be able to understand in more detail the behavior of the critical 
properties as a function of d. Eq. (p9D, to the best of our knowledge, does not lend itself 
to an exact treatment for arbitrary values of /3 and d. However, on the critical surface for 
weak-strong transition, we find that it is possible to turn Eqs. (|3^ ) -(l40| ) into a homogeneous 
(eigenvalue) equation. For c? < 4, this eigenvalue problem can be solved analytically. For 
c? > 4, the problem can be solved numerically with great precision. 



It is also worth pointing out that at criticality for (i < 4, Eq. ( p9D can be solved by 
applying a method of Gaudin and Kostov for the study of 0[n) spins on random surfaces. 
There is in fact an exact mapping of our weak-coupling critical saddle point equation to 
that of Ref. 



15 



with (i = n + 2. However, for c? > 4, the Gaudin-Kostov's solution become 
pathological. In contrast, for the simplicial models, we find that a consistent solution exists 
for ci > 4. 



The determination of /^c can be achieved by considering Eq. ( |39D in the limit when the 
equality in the weak coupling constraint, Eq. (^Tj), is reached. On the weak coupling side of 
criticality the condition p(0) = implies 
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arbr. 



1 d-2 f^- dy 



pc{y) 



2 J-c y J(^bc + y){y + ac) 



111) 



There are two possible solutions, (i) ac = 0, and (ii) ^ 0, with 



y J{bc + y){y + a,) 



:ii2) 



If Eq. (^) is solved by = at /3 = Eqs. (|3|) and (|0D both reduce to 

1 d-2 fbc dy p,{y) 



Pc{z) 



{be - z)z 



TT 



(113) 



This indeed applies for < < 4, and we shall find an explicit solution of Eq. ( p.l3|) , which 
agrees with a result previously found by Gaudin and Kostov |T^. In particular, we find that 
(3^ = 1/d for < ci < 4. 



When d > 4, Eq. (|111|) must be solved with a^, > 0. Substituting Eq. (|112|) into 
Eqs. (p9|)-(|i0|), one arrives at a homogeneous equation 

r^c dy pe{y) 



Xpc{z) = {be - z){z - ac) 



»^ y(y + ^) ^{bc + y){y + ac) 



(114) 



where A = Eq. ([Tl4D can be solved numerically with great accuracy. Surprisingly, 

the relationship (3c = l/d was found to be satisfied within machine precision. The essential 
feature of the solution for > 4 is the strong coupling relationship = ^ 0, while in 
weak coupling necessarily r = 0. Since Pc{z) is the same on both sides of the transition, at 
criticality one finds from Eq. (^) 

^ ' , (115) 



1 P.i^y-c 



4/32 7a. d{d - 1) A{d - 1) 



and as a consequence a first order phase transition is observed. 



A. Critical Solution for < d < 4 



When < d < 4, Eq. (^9]) at criticality can be solved by assuming that a^. = 0, 
as suggested by our analytic results discussed in Section III. Let us therefore focus on 
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Eq. (|113|) . This equation on the first sight suggests that p{z) would vanish at 2; = as y/z. 
However, it is easy to verify that, upon substituting this behavior into the right hand side 
of the equation, this square-root behavior is in fact inconsistent. Based on our earher exact 
analytic solutions, we assume that Pc{z) vanishes at z = faster the ■\/z; it follows that the 
square-bracket in Eq. ( p.l3| ) must also vanish at z = 0. As a consequence, we have 



and we again arrive at a homogeneous equation 

AM.) =.»'^V^r^T7l#=. (117) 

^ Jo y + zy■^|^^Jbc + y 



Note that this homogeneous equation connects smoothly with that appropriate for ci > 4, 
Eq. (ini), with = 0. 



It is convenient to change variable from z to uj, uj + uj ^ = 2hc/ z. In solving for uj in 
terms of 2;, we shall choose the branch 1 < < 00 so that Eq. (|1 1 7|) becomes 



Jl [uj' + Uj)[ljj' + UJ ^) 



where f{uj) = {1/2){lo + u-^)p,{z). 

Although f{uj) is originally defined only for the interval 1 < < 00, the right-hand 
side of Eq. ( p.l8| ) provides a natural extension to the region Q < uo <1. With this extension, 
one finds that, over the positive axis, < c<j < 00, 

/(.;) = -f{oj-^) . (119) 



It is then straightforward to verify that this extension can also be made for the right-hand 
side of Eq. ( |118| ) so that it becomes 

^/ ^ = ^ / I ,\ ,^ Jj i^) ■ 120 

2 Jo [uj' + uj)[uj' + uj 



Let us next treat Eq. (|120D as an eigenvalue problem, and consider the Ansatz where 
/(a;) = c[lj^ — uj~^]. Using the technique of contour-integration, it is easy to verify that 
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this indeed is an eigenvector with eigenvalue A = vr/ cosirO. Since A = it follows that 
d = 4cos^^. With < d < A, one has 9 real and < < 1, which allows a solution where 
Pc{z) is positive definite! 



Using the normalization condition for together with the criticality condition 

2 



/q^" dz' / z' Pc{z') = 1, we can fix the normalization constant c = ^cos^ and the end point 



be = ftg^ . One then obtains 



, , cos^sh6'M 



where = uj and < m < oo. One can show that Eq. ( |121|) reproduces the known critical 
solution when d = 2 and 4. When substituted into the critical equation at d = 3 Eq. (|121|) 
is numerically found to be a satisfactory solution. 



To determine the critical value /5c, we can re-express Eq. ( |116D in terms of f{uj) as 



^ ^ ' rfo; ^ ^ , , V H- (122) 



(3c V2bc Jo 1 + ^ 



Again, by an contour integration, one arrives at the remarkable result 



B. Criticality for d > i 



The solution discussed above does not apply to the case d > A because the analytic 
continuation of Eq. ( |121| ) for the critical density would no longer be positive-definite in the 
interval [0, be]; we must choose the alternative, Oc 7^ 0. We have previously seen, with Oc 7^ 0, 
how the criticality condition for Pc, Eq ( |112|) , and the homogeneous integral equation for pc, 



Eq. ( p. 14| ) , can be obtained, approaching from the weak coupling regime, by enforcing the 
condition p(0) = with r = 0. It is instructive to see how these equations can be similarly 
derived from the strong coupling regime. 

Starting with Eq. (|40|), one is working within the strong coupling regime where the 
constraint, Eq. (pOD, is automatically satisfied with a = a/t 7^ 0. It can be shown that the 
criticality condition, Eq. ( p.l2| ), corresponds to a situation where a zero of p{z) enters at 
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z = a. That is, as one increases j3 beyond j3c, the positivity of p{z) would be violated, thus 
terminating the validity of the strong coupling solution. As pointed earlier, with Eq. (|112D , 
the strong coupling equation, Eq. (^OD, again leads to Eq. ( |114| ). As a consequence, for d > 4, 
when one approaches (3c from the strong coupling regime, one finds that Oc = ^/rc 7^ and a 
first-order phase transition occurs. The solution to Eq. ( [L14| ), on the critical point, subject 
to Eq. (pOD and Eq. (^2]), can be found numerically. In order to have a better behaved kernel 
when d is close to 4 we define a new function H, 

Hc{z) = P^^'^ . (124) 

zJ {b — z)[z — a) 



With A = the integral equation Eq. (114), and the constraints, Eqs. 
equation which determines /3c, Eq. (|112|) , become 



XHJz] 



^ dy ^{h-y){y-a) 

j^c\y) I 

-y + z J{h + y){y + a) 



2|), and the 



(125) 



f dz Hc{z) zJb - y){y - a) = 1 
[ dz Hc{z) Jb-y){y-a) = 1 , 

J a 



(126) 
(127) 



1 r'' a/(6 - z)(z- a) 

- = 2{d-2) dzHciz) 



^{b + z){z + a) 



(128) 



The solution of the integral equation Eq. (|125|) can be done numerically by discretizing 



the kernel. After the discretization, the problem is reduced to an eigenvalue problem of a 
real non symmetric matrix. There are several ways to discretize the kernel. Any rule of 
numerical integration is a discretization rule for the kernel. It is known that for integral 
equations the best discretization rules are the Gauss quadrature rules [^]. There are several 
Gauss quadrature rules. We used the simplest possible: the Gauss-Chebyshev rule. All these 
rules require to map the integration interval to [—1, 1]. Thus we perform the following change 
of variables 

b — a a + b , ^ 

^ = + ^ , (129) 



or 
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2 a + b 

C = 7 ^ - 7 ■ 130 

— a — a 



A^c(c) = I . , 7, I:...: , (131) 



Under this change of variables the integral equation becomes 

dC ^'^^^ ^^^^ 
' ^ + C + 2l?^ v/(e + f^^)(e + ^) 

where k = a/b. The constraints take the following form 

J\ di H,(o yn; = ij^?^ , (133) 

and the equation for Pc is 

1 ./I - £^ 

— = 6(ci-2)(l-/s:)/ ^ ^ ^ . (134) 



The solution to Eq. ( |131| ) can be found up to an overall constant C, (assuming that A 
is a nondegenerate eigenvalue). This constant and the upper bound, b, of the support of Pc 
can be computed using the constrains Eq. ( |132| ), Eq. ( |133| ). Let's denote an eigenfunction 
of ( |131| ) by He- Then He, which is the function that is positive in [-1,1] and satisfies the 
constraints, ( |132| ) and ( |133| ), is related to He by H^ = CHc- If we now define 



h = I d^ HeiO + - ' (135) 

1 L Ki 



d^HeiO^l-e, (136) 
-1 ^ 



then 



u ^2 2 
b = — 



Jl 1 - K ' 



r - ^ 

-'2 



(137) 
(138) 
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From the above formulas it is obvious that one can fix k, solve Eq. (|131|) and then find 
the eigenvalue A which has a positive definite eigenfunction. Because the problem is well 
defined one expects that there exists only one such function. This expectation is confirmed 
by the numerical results. It turns out that the eigenfunction with the largest eigenvalue is 
the one which is positive definite in [—1, 1]. The nth eigenfunction has n — 1 zeros in [—1, 1]. 
Thus for a given n one computes b{K,) , a{K,) , A(/t) , Pdn) , d{K) = 27iX{k)^^ + 2. 



Using the above numerical method we have computed ac, be with great precision for 
d in the interval (4.4,250). Combining the numerical results with the analytical for d < 4, 
and be are plotted in Fig. as functions oi 1/d for d > 2. The functions add), bc{d) are 
continuous functions of ci at ci = 4. 



Several interesting features now emerge from an analysis of this data. On the one hand 
we can fit the functions add), bdd) with great accuracy as power series of d — A around d = 4 
and they agree with the corresponding weak coupling expressions up to very high orders. 
On the other hand, if one does a careful extrapolation of a^, be to d = 4 a new feature is 
seen ( Figure. |^). The upper limit, bdd), extrapolates linearly in (i — 4 to vr, consistent with 
analyticity in the d — 4 series expansion. The data alone determines the intercept &c(4) to be 
IT to an accuracy of 10~^. When one examines the lower limit ac, it also approaches zero as 
d approaches 4. However it does not go to zero as a simple power. The more we improved 
our data near d = A the higher the effective power became. It appears that ac may have 
an essential singularity at c? = 4 vanishing faster than any power. Since the discontinuity 
of the internal energy on the first order line is given by da1/{4:{d — 1)), this is pertinent to 
the critical properties at the end of the first order transition. The log-log plots of ac{d) and 
vr — bc{d) in Fig. |^ clearly support these observations. 



In Fig. 1^ we have also been helped by an expansion of the functions (oc + 6c)/2 and 
ajjc in powers oil/d. The coefficients of this expansion have been determined by best fits on 
the numerical results and found to be consistent with integer numbers within the precision 
of our determination. This result is also consistent with the results of the 1/d expansion 
which we shall discuss in the next section. In particular we found 

ac + bc _ 2 2 _ 6 _ 20 _ 48 _ 92 _ 118 _ 42 

~^2~~ ^ d^ d^~ d^~ d^ ~ d^ ~ d^ ~ 1^ ~ d^^ "" ^ ^ 

12 40 8 40 128 328 694 1112 , , 



The last integer term in both equations is uncertain. 



Furthermore from the numerical computation of the f3c for d > 4, we can see the d < 4 
result /?c = ^/d still holds above the critical point. We performed the numerical calculations 
in double precision and we see no deviation at all from the 1/d law. The deviation of 
|/3c — l/d\ from zero is determined to be less than 10^^^ for d in the interval (4.4,250). 
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V. SUMMARY AND CONCLUSION 



The 1 /N expansion of matrix models has recently been used as a discrete representation 
for summing over random surfaces and, through the "double-scaling" limit, for studying low- 
dimensional string theories. Even more importantly, the large-A^ expansion has provided 
us with a scheme for addressing non-perturbative issues in non-Abelian gauge theories. For 
instance, many qualitative features of QCD, e.g., confinement, the OZI rule, etc, can best 
be understood in a large- iV setting. 



However, quantitative progress in these directions has been slow due partly to the 
technical difficulties associated with the large number of independent "loop" variables in 
this limit. Nevertheless, it has been possible to gain useful insights into various interesting 
situations by utilizing as guides solvable models involving a small number of matrices, e.g., 
models involving two Hermitian matrices. Much less is known for models involving unitary 
matrices. Our current work not only adds to the list of solvable models in this category but 
also introduces new techniques for addressing matrix model studies in the large- iV limit. 



In this paper, we have studied the large- structure of simplicial chiral models defined 
on a (i — 1 dimensional simplex as one varies d and the coupling (3. By exploring the global 
U (N) X U (N) symmetry and by introducing an auxiliary complex matrix field, we are able 
to reduce the problem to that of solving for the eigenvalue distribution of a single Hermitian 
semi-positive-definite matrix in the large- limit. In addition to providing exact large- 
solutions for several specific values of d, we are able to identify and solve the strong-weak 
criticality problem for all values of d, < d < oo. 



For < d < A, analytic solutions for pc can be found. Interestingly we find that 
the criticality occurs precisely at /5c = 1/^, as suggested by our previous studies [|TU|. We 



find that the transition is third order. For 4 < d < oo, the criticality can also be studied 
by solving a homogeneous integral equation. However, we are only able to carry this out 
numerically. Within numerical accuracy, we have shown that the criticality again takes place 
ai = ^/d, but with a first order transition. 



Since we are able to reduce a d — 1 dimensional simplicial chiral model to a model 
involving a single complex matrix, with d entering as a parameter in the effective action, the 
large- limit can thus be solved by finding a density function for the eigenvalue distribution. 
Unlike usual matrix models where all eigenvalues lie in a single connected band in the large- 
A^ limit, this model effectively involves two bands, a "right-band" where > and a 
"left-band" where Xj < 0. However, unlike other two-band problems p9[, the distribution 
over these two bands are correlated. This new feature presents a challenge which cannot be 
handled by a conventional large- A^ treatment. 
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Our key result is the reduction of the above problem to that of solving a single inhomo- 
geneous integral equation for the eigenvalue distribution of a single Hermitian semi-positive 
definite matrix. Although we could not find a closed form solution to this equation for arbi- 
trary d, we are able to solve it in several interesting special cases and we set up a systematic 
numerical approach to the solutions. We have found that the critical surface is defined by 
I3c = 1/rffor all d. 

For small d, < li < 4, the models exhibit the third order Gross Witten transition. 
Indeed for c? = 1, 2, 3 they coincide exactly with the chiral chains studied earlier by Brower, 
Rossi and Tan 0. In this region, the criticality is related to that of 0{n) spin models on 
random surfaces, as discussed by Gaudin-Kostov |]15[. For > 4, however, there is a first 



order transition ending at the "upper critical" dimensions d = A. It therefore appears that, 
from the perspective of the double- scaling limit, the most interesting situation corresponds 
to < d < 4. We have found that the point d = 4, having a logarithmic singularity, 
corresponds to a = 0. In the language of the double- scaling limit, this corresponds to 
having a vanishing "string susceptibility", a = 'jstring = 0, where Csing ~ (/5 — l3c)~^"*""^, 
which formally correspond to that resulted from a c = 1 CFT theory. This calls for further 
studies around d = 4 which can provide further insight into possible different mechanisms 
for generating c = 1 physics. One way is to vary d near d = 4. Another approach is to stay 
ai d = 4, and embellish the model by relaxing the "permutation symmetry" of the original 
d = 4 simplicial chiral model. This will be presented in a subsequent publication. 



ACKNOWLEDGMENTS 



We are deeply indebted with Prof. G. Cicuta for bringing Ref. |T5[ to our attention 
and for useful conversations. This work was supported in part by the U. S. Department of 
Energy, under grant DE-FG02-91ER400688, Task A. 



APPENDIX A: CRITICAL PROPERTIES OF CHIRAL CHAINS FOR L = 2, 3, 4, oo 



Chiral chain models are defined by the partition function 



i=l 



'NpY.Mu^ul, + uJu,^^] 



i=l 



(Al) 



with Ui = Ui+L- Free energy density, internal energy and specific heat are given by 
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Fl 
Ul 
Cl 



1 



2 d(3 



In Zr 



0- 



d/3 



(A2) 



When L ^ oo, Zl can be reduced to the partition function of the Gross- Witten 
single-hnk problem ||^ 



dU exp \N(3Tt{U + U^) 



(A3) 



thus sharing the same thermodynamic properties. The free energy density at = oo, 
F = ^InZ, is piecewise analytic with a third order transition at /5 = /?c = 1/2 between 
the strong coupling and weak coupling domains. The large- A^ limit of the specific heat is 

Coo = P^ for P<(3,, 

Coo = \ for (3>(3,. (A4) 



The behavior of Coo around (3c can be characterized by a specific heat critical exponent 
a = —1. It is worth noting that an analysis of the double scaling limit, N ^ oo and (3 ^ f3c, 
allows the determination of the correlation length critical exponent, z/ = 3/2 |l7| , p!8[| , and 



that a and u satisfy a hyperscaling relationship associated to a two-dimensional critical 
phenomenon, 2i^ = 2 — a. This fact is related to the equivalence of the double scaling limit 
with the continuum limit of a two-dimensional gravity model with central charge c = —2. 



It has been shown that in the context of single-matrix models the parameter A^ plays 
a role quite analogous to the volume in ordinary systems, and double scaling limit turns out 



to be very similar to finite size scaling in a two-dimensional critical phenomenon ||19| , p0 
As a manifestation of this fact, it has been observed that in the Gross- Witten single-link 
problem, (i) the asymptotic approach of the complex Z{N,P) zeroes closest to /?c, P{N), 



toward the real axis occurs at a rate determined by the correlation length exponent |21|] , 

lmp{N) oc A^-^ , (A5) 

and, (ii) for sufficient large A^ the position of the peak of the specific heat, Ppeak{N), behaves 
as 

0eak{N)-(3c OC N-^ . (A6) 
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We recall that in ordinary critical behaviors finite size scaling leads to relations of the type 
(^-A6) with replaced by the size of the system. 



The L = 2 chiral chain again corresponds to a Gross- Witten model with (3 replaced 
by 2/5, thus obtaining jSc = 1/4 and the same critical exponents. 

Solutions of the models with L = 3,4 have been found in Ref. J^. The results we 
present in the following without details on the derivations are easily obtainable from the 
analysis of Ref. . 

We recall that the chiral chain with L = 3 is equivalent to the simplicial model with 
d = 3. In this case f3c = 1/3 and the phase transition is still third order. In the weak 
coupling region, P > /3c, the N = oo specific heat is given by 

Therefore close to /3c, 

Ct^ = i - ^(/^ - ^^)'^' + 0{P-Pc). (A8) 

Similarly in the strong coupling region, /? < /?c, and close to /?c, 

Ct^ = ^ - ^(/5c - f^y/' + - (3) . (A9) 



Then the strong and weak coupling expressions of C3 show that the critical point /5c = 1/3 
is third order and a = —1/2. 



For L = 4 the study of the critical behavior around jSc = vr/8 is slightly subtler. In 
the weak coupling domain the N = 00 internal energy can be expressed as 



Uf^ = 2/? - - - 2/35^ 



1 



(AlO) 



where 6 is implicitly determined by the equation 

^ [e (VT^) - 6'F (Vl^ 



1 . 



(All) 
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(This equation comes from the normahzation of the eigenvalue distribution p{9) introduced 
in Ref. 0). Since 5 = at /?c = ^r/S, in order to study the critical behavior close to j3c we 
expand Eq. ( |A11| ) around 5 = 0, obtaining the following relation 



and therefore 5"^ ^ (3 — j3c apart from logarithms. Furthermore we have 

d5^ 16 



d(3 7rln(4/5) 
We then obtain for the specific heat 

TT^ 1 



+ 0{6') . (A13) 



when /? — > 

A similar analysis can be performed in strong coupling, where the internal energy can 
be written as 

Ut^ = 2/5 - ^ + 2/5^ , (A15) 



with ( implicitly defined by the equation 

8(3 1 



E Jl-e] = I. (A16) 



Since C = at jSc, we expand Eq. ( [A16| ) around C = obtainin 



- \'(^4 + l]+0(C) . (A17) 



Pc 2^ V C 2 

Consequently, ~ /?c — apart from logarithms and 
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We then obtain when (3^/3^, 



Cr = — + , , , , + 0((3c - (3) . (A19) 

^ 32 8 161n(4/C) 



A comparison of Eqs. ( [A14| ) and ([A 191) leads to the conclusion that the phase transition 
is again third order with a critical exponent a = 0~ . The critical exponent u could then be 
determined by using the 2-d hyperscaling relationship, obtaining u = 1. This value of v was 
confirmed by a numerical Monte Carlo study of the scaling of the specific heat peak position 
at finite A^, we indeed observed a behavior like Eq. ( |A6|) compatible with u = 1 within a 
few per cent of uncertainty. 



In Fig. ^ we plot the specific heat versus f3 for L = 2, 3, 4, oo. In conclusion we have 
seen that L = 2, 3, 4, cxd have a third order phase transition at increasing critical values 
Pc = i) |) f ! with specific heat critical exponents a = — 1, — 1/2, 0~, — 1, respectively. 
Notice the behavior of a with respect to L, which is increasing for L = 2,3,4 reaching the 
limit of a third order critical behavior, but then in large-L limit it returns to a = — 1. 



APPENDIX B: STRONG COUPLING EXPANSION OF CHIRAL CHAIN 

MODELS 



Strong coupling series of the free energy density of chiral chain models are best gen- 
erated by means of the character expansion, which leads to the following result 

Fl{/3) = F{(3) + , (Bl) 



where F{f3) is the free energy of the single unitary matrix model {F{j3) = In Z and Z is 
given by Eq. ([A^)), and 

^ ln^4)^(,)(/5)^ (B2) 



(r) 



X](r) denotes the sum over all irreducible representations of U{N), d(^r) and -2(r)(/5) are the 
corresponding dimensions and character coefficients. The calculation of the strong coupling 
series of Fi{l3) is much simplified in the large- A^ limit, due to the following relationships 

F{(3) = 13^ + , (B3) 
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and 

= + , (B4) 



where is independent of P and n is the order of the representation (r). Exphcit expres- 
sions of d(^r) and Z(^r) are given in Ref . |^ . Notice that the large- strong couphng expansion 
of Fl{(3) is actually a series in {3^ , 



Y.c{n,L)(3-'^ . (B5) 



Fl represents also the generating functional for the "potentials" W{n, L) introduced in 
Ref. , in the context of the strong coupling expansion of more general models, indeed the 
following relationship holds 

Win,L) = |c(n,L). (B6) 



It is important to recall that the large- character coefficients have jumps and singularities 
at /5 = I , and therefore the relevant region for a strong-coupling character expansion is 
/3<i. 

We have analyzed the strong coupling series of chiral chain models in order to inves- 
tigate their large- iV critical behaviors for L > 4. Given the simple behavior of the large-iV 
limit of F{j3), we considered only the contributions from F{(3), thus working with series 
in (3^ . We generated about 15 terms for each L < 10 and analyzed, as series in the 
specific heat derivative, which diverges at the critical point in a third phase transition. We 
employed the integral approximant technique |2^,24,2^, which at present seems to be one 



of the most powerful method of resummation. In particular we considered integral approxi- 
mants obtained from first order linear differential equations. 

Let us begin with the results obtained for the known cases L = 3, 4. For L = 3 already 
15 terms in the series (in (5^) suffice to get (3c = \ and a = —\ with a precision of about 
10~^ and 10"'' respectively. However it is worth noticing that in the analysis of the specific 
heat derivative we found spurious non-diverging singularities on the positive real axis for 
P<(3c. 

Concerning the L = 4 case, it is known that the integral approximant resummation 



analysis cannot reproduce an a = 0~ singularity type [0] and therefore it is not really 
suitable to this case. Anyway, we obtained a good determination of /?c, we found | up to 
about 10~^, and a rather stable but wrong exponent, a ~ —0.18, which should somehow 
simulate the logarithmic corrections found in the Appendix A, given that they cannot be 
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generated by the differential equation solution. Again we found spurious non-diverging 
singularities for (3 < P^- 



The strong-coupling analysis starts giving new information when L > A. Due to the 
persistent presence of spurious singularities, guided by the L = 3,4 analysis, in all cases we 
considered the first diverging singularity on the positive real axis as an estimate of the true 
critical point. For L = 5 we obtained quite stable results: Pc — 0.43756 and a ~ —0.17. We 
should say that the L = 4 analysis suggests some caution in accepting this estimate of a, it 
could still be a masked a = 0^ . The analysis of L = 6 series gave a rather stable estimate of 
the critical point f3c — 0.4737, but unstable exponents (although negative and small). Similar 
results were found for L > 7: ~ 0.504 for L = 7, ~ 0.526 for L = 8, ~ 0.57 for L = 10. 
Notice that, unlike the L < Q cases, these values cannot be considered as an estimate of the 
critical point. They are indeed larger than |, that is out of the region where a strong coupling 
analysis can be predictive, and therefore something else must happen before, breaking the 
validity of the strong coupling expansion. An example of this phenomenon comes from 
the Gross- Witten single-link model (recovered when L — > oo), where the strong coupling 
expansion of the N = oo free energy leads to an analytical function not having singularity 
at all, F{P) = thus Pc = \ cannot be determined from a strong coupling analysis. 

Of course we cannot consider this analysis satisfactory, but from it we may hint at a 
possible scenario. As for L < 4, for L = 5, 6, that is when the estimate of Pc coming from 
the above strong coupling analysis is smaller than | and therefore acceptable, the term F{P) 
in Eq. ( PT| ) should be the one relevant for the critical properties, determining the critical 
points and giving a ^ —1 (maybe a = 0~ as in the L = 4 case). For L > 7 the critical 
point may not be a singular point in strong or weak coupling, but just the point where 
weak coupling and strong coupling curves meet each other. This would cause a softer phase 
transition with a = — 1, as for the Gross- Witten single-link problem. We expect Pc < \ also 
for L > 7. 



This scenario would be consistent with the analysis of Green and Samuel |]2^, who 
studied the behavior of the link determinant (i.e. {deiUiUlj^i)) to determine the critical 
points. The values of Pc we found for L = 5, 6 are consistent with their estimates. 



APPENDIX C: d = 4 MODEL: WEAK AND STRONG COUPLING EXPANSION 

AND SERIES ANALYSIS. 



We have briefly discussed in Section |1I1B| the possibility of performing weak and strong 
coupling expansion in the d = A model starting from an expansion in the powers of k"^. Here 
we want to give more details on the concrete implementation of this program. 
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Let us first of all for convenience define a few auxiliary functions of k"^, the labels s 
and w are to remind the strong and weak coupling expansions, which we shall treat on the 
same footing. In terms of standard elliptic integrals we define 



n„ik') = - / ^7^=^dC = -E{k) , (CI) 



n3(P) = - / ^^T^^rfC = - \E{k) - (1 - P)ir(A;) , (C2) 



vr Jo a/1 — vr 



dw(A;') = -K{k)n,{k^) , (C3) 

TT 



= -K{k)k^n^{k^) . (C4) 

TT 



let us introduce the function 

•2\2 rk 



Pik') = (-) v'F-C2y^^^(C^A;K, (C5) 



and notice that, because of the properties of elliptic integrals and of Eq. ([7^) we may obtain 
from Eqs. (|70l-^ 

- d{k^)-pm ^^^^ 



holding for the proper choice of indices both in the weak and in the strong coupling phase. 
In order to set up our expansion we therefore need power series representation for the 
functions ^E{k), ^K{k) and piji^). The elliptic integrals of the first and second kind have 
simple known expansions. 



2 



oo 

K{k) = 1 + ^co„(A;2)" , (C7) 



^ n=l 



where 
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_ (2n!)^ 



The elliptic integral of the third kind, needed in the construction of admits the 

following expansion in the powers of the first argument, 

r\ OO 

-^(C^A;) = EC^'^K^'), (CIO) 



where 



/^A;') = - / ^ (Cll) 

Trio 

satisfy the recursion equation 

{21 + - 21(1 + e)Ii{k^) + {21 - 1)/,_i(A;2) = , (C12) 



with boundary condition, 

h{e) = -K{k) , 

TT 

h{k') = ^ [i^(fc) - E{k)\ . (C13) 



In turn we may show that 

Mk') ^ f'^^^'^/l^'CdC (C14) 

are related to by 

Ji{k') = Ii{k') - (1 + k')Ii+,{k') + k'li+2{k') . (C15) 



As a consequence the function p{k^) is completely determined once Ii{k^) are known, by the 
relationship 



00 



p{k') = Y.^iik')Ji{k'){ky+' . (CI6) 

1=0 
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The recursion Eq. ( |C12|) may be solved if we expand the functions Ii{k'^) in a power series 
of k\ 

oo 

Iiik') ^ J2^inik'r. (C17) 

n=0 



It is tedious but straightforward to verify that 

[(n + /)!n!22"+f ^ ^ 



satisfy the recursion. 

One may also for convenience define 

oo 

Jiik^) = Y^dUk'r, (C19) 

n=0 



and finds from Eq. ( |(J15| ) that 

~2(2n-l)(n + / + l 



din = -7777^ , , , • (C20) 



Direct substitution of the series expansions thus obtained in Eq. ( PBD allows to construct an 
expansion of (3 (or l/j3 respectively) in powers of k"^ which a simple symbolic manipulation 
program can easily extend to extremely high orders. 

In order to make our analysis complete we must manage to extend our discussion to 
the evaluation of a physical observable. By recalling Eq. (0), we notice that the expression 
for the internal energy in the d = 4 case is 

U = t dz' p(z') iz'^ -r)---\, (C21) 

48/52 Ja ^ 3 12/5 ' ^ ^ 



which gives rise in weak coupling to 



while in strong coupling we obtain 
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3ns(P)2 Jo ^^^^^^^ ^^33 ns(P) 



By considering the explicit form of the functions Pw(C) ^'iid Ps(C) which we obtain from 
Eqs. ([70| - |7TD , it is straightforward to parameterize U by 



^ 4 fcV(fc2) _ 1 _ 2 rf(P)-p(P) 

3 [rf(A;2) -p(A;2)]n(A;2)2 3 3 n(A;2)2 ' ^ ^ 



where the functions r(A;2) can in turn be reconstructed, by using the results presented in 
this Appendix, in terms of the functions Ii and Ji. The results for the two regimes are 



(e) = ^ Jiik') - {ey (C25) 



1=1 



and 



rsik') = Y: Mk') [Ii-i{k') - k'liik')\ (C26) 
1=1 



respectively. By expanding Eq. ( P24|) in power series of k'^, inverting Eq. (|C6D and substi- 
tuting A;2 as a function of or /3 respectively, we obtain the standard weak and strong 
coupling series for the internal energy. 

In conclusion we present some terms of the strong and weak coupling series of the 
internal energy, obtained by implementing the procedure outlined in this appendix: 

U,= (3 + 2(3^ + 2(3^ + 4(3^ + 28 (3^^ + 38(3'^ + 8 + 440 f3^ + 1936 /3^° + 1712 f3^^ 
+4160 + 62160 + 178072 + 101038 + 1215704 P^^ + 9259720 P^'^ 
+17052880 p^^ + 15519376 p^^ + 291277184 p^^ + 1351546592 p^^ + ... (C27) 



AP [Apy^ (4/ij^ (4/3j* [4:P)° [AP)'' 

nnniQ7i nnniARQ nnnii5^5^ n nnnsQS n nnn797 n nnnfinn 

+ ... (C28) 



The strong coupling series was also generated by using the more general approach of Ref . 
finding the same results. 
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APPENDIX D: NUMERICAL APPROACH BASED ON LARGE d EXPANSION 



The Ansatz (^) can also be used as a starting point for numerical approximations 
based on the very simple consideration that real analytic functions of ( can be approximated 
with any assigned precision by polynomials of sufficiently high degree in ( itself. We may 
therefore introduce the n-th truncations of Q{C) and R{C) respectively by the definitions 



Qn(C) = E^^C, 

i=l 
n 

RniC) = J2^^C■ 



(Dl) 



i=l 



We may now explicitly perform a Laurent series expansion around the point ( 
form, 



in the 



dQniC) 



An) 
•-1 



(n) 
C2 



(D2) 



i=l 



where the coefficients pi, po and are completely determined in terms of qi, ci and C2. 



We may now notice that the asymptotic condition on 0n(C) forces us to impose the 
relationships 



Pi + Vi = i > 
Po = , 
P-i = I . 



(D3) 



At this stage i?„,Ci"\c2"^ are completely determined in terms of the coefficients qi and the 
parameter z^^\ We may now consider the effect of substituting (pniC) i^ito Eq. (|96D and 
power-series expanding in if we define (l)l^\—2z) to be the k-th derivative of the func- 
tion (pn evaluated at the point —2z, we may turn Eq. (^6]) into the following approximate 
relationship, 



= 2[i?„(C) + rfQ„(C)] + (rf-2)f:^^0i'=)(-25) 



k=0 



k\ 



(D4) 



which in turn decomposes into n + 1 equations in the n + 1 unknowns g^, 
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z 



(d-2)0„(-25) 



2/3 



d 



1 




2/9 







\ <k <n . 



(D5) 



These equations may be solved numerically for arbitrary values of /3 and d and offer bet- 
ter and better approximations to the true eigenvalue distribution (by applying Eq. ( |109|) 
with increasing values of n). Numerical experiments have shown that when (5 '^^Xjd even 
extremely small values of n give quite accurate predictions, while around the criticality, 
corresponding in this language to the condition 



which determines (5c-, but the accuracy is definitely weaker. This fact did not prevent us 
from determining the location of criticality, by the use of n < 8, in the cases d = 3,4,5, 
with a precision of about 1%. 

In the context of this discussion it is important to observe that (even approximate) 
knowledge of 0(C) implies an (approximate) knowledge of the moments of the eigenvalue 
distribution, which may be obtained from the Laurent series expansion via the relationship 




(D6) 




OO -1 „^ OO 



(D7) 



which in particular implies that 




(D8) 



and in turn we may extract the internal energy via the relationship 



d[d -1)U = — [z^ + 2zp.2 + P-s) -d- 



(D9) 
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FIGURES 

FIG. 1. dC versus d/3 for the d = 2, 3, 4, oo simplicial models. 
FIG. 2. Oc and be versus 1/d. The crosses mark the points corresponding to d = 4. 
FIG. 3. Log-log plot of Uc and it — he versus d — A. 
FIG. 4. Specific heat versus (3 for the L = 2, 3, 4, oo chiral chain models. 
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